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Abstract 

We propose a mechanism to explain the fluctuations of the ground state 
energy in quantum dots in the Coulomb blockade regime. Employing the ran- 
dom matrix theory we show that shape deformations may change the adjacent 
peak spacing distribution from Wigner-Dyson to nearly Gaussian even in the 
absence of strong charging energy fluctuations. We flnd that this distribution 
is solely determined by the average number of anti-crossings between consec- 
utive conductance peaks and the presence or absence of a magnetic field. Our 
mechanism is tested in a dynamical model whose underlying classical dynam- 
ics is chaotic. Our results are in good agreement with recent experiments and 
apply to quantum dots with spin resolved or spin degenerate states. 

PACS numbers: 73.23.Hk, 05.45.-Fb 
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Most phenomena in mesoscopic electronic transport are well understood invoking the 
picture of a single quantum particle moving coherently in a disordered or complex potential. 
In the case of semiconductor quantum dots, the complexity of the background or confining 
potential is manifest in the chaotic nature of the underlying electronic classical dynamics. A 
large amount of literature indicates that the statistical properties of the spectra of chaotic 
systems are successfully modeled by random matrix theory (RMT) 0. When applied to 
transport in open ballistic quantum dots of irregular shape, this theory predicts a strong 
universal behavior for the statistical fluctuations of the conductance, similar to that found 
in disordered systems |Q. Over the last years, several experiments confirmed this behavior, 
making open electronic cavities one of the paradigms of quantum chaos @,|^. 

Recent experiments showed that RMT is also very successful in predicting the 
distribution of Coulomb blockade peaks heights in nearly isolated, many-electron quantum 
dots at low temperatures and zero bias P]. In this case, the theoretical description is based 
on the so-called 'constant interaction' (CI) model 0, where the electron-electron interaction 
is taken into account through a fixed (capacitive) charging energy term in the Hamiltonian. 
The electronic many-body wave function is then formed by filling single-particle eigenstates 
whose fluctuations reflect the chaotic nature of the classical motion within the dot. 

However, new data related to the many-body ground state energy fluctuations in these 
systems are at variance with the standard RMT predictions. It was observed that adjacent 
peak spacing energies do not fluctuate according to the celebrated Wigner surmise |[T0HT2|. 
Moreover, the existence of strong correlation among the heights of nearby peaks and the 
apparent lack of spin degeneracy seemed to indicate that other effects, such as wave function 
scarring |13[, interaction- induced correlations, and charging energy fluctuations |]T0|,O,ll5 



also play important roles. 

The aim of this letter is to show that, although effects beyond the single-particle ap- 
proach may be relevant to the general understanding of quantum dots, there is yet a simple 
mechanism based on the CI model and RMT which explains the basic feature found in 
Refs. ||10Hl2|, namely, the approximately Gaussian form of the peak spacing distribution. 
The essential point which has been missing in all previous analyses of the problem is that 
the gate voltage, swept to produce a series of Coulomb blockade peaks, also continuously 
deforms the confining potential seen by the electrons in the quantum dot [|l^ . As a conse- 
quence, ground state energies, which are measured from the positions of consecutive peaks, 
correspond to different parametric realizations of the one-body Hamiltonian. Since adjacent 
energy eigenvalues tend to lose correlation rapidly as their parametric distance increases 



r?| , the experiments in question are likely to be measuring the spacing distribution of un- 
correlated chaotic eigenstates. If the dot shape dependence on the gate voltage is sufficiently 
strong, we find that the peak spacing distribution approaches a Gaussian, regardless of the 
existence or not of spin degeneracy in the single-particle spectrum. 

In order to present the quantitative aspects of our argument, we begin indicating the 



condition for the occurrence of a conductance peak within the CI model [18 



En + nU- er]Vg = Ep, (1) 

where En is the n*^ single-particle energy level (n ^ 1), ?7 is the dot charging energy, Vg 
is the gate potential, and Ep is the Fermi energy (chemical potential) in the leads. The 
coefficient 77 is a function of the capacitance matrix elements of the dot 
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In all experimental setups where plunger gates are used to define both shape and depth 
of the confining potential, one expects that variations in the potential Vg will continuously 
deform the dot. This is the actual situation at least for the experiments in Refs. [|11|,0. To 
establish the direct relation between Vg and the dot shape we need to know the experimental 
set up in detail. While such study is desirable, it is not crucial for our semi-quantitative 
analysis and will not be pursued here. We therefore simplify the problem parametrizing the 
shape deformations by a generic variable X, such that En = En{X). 

In view of Eq. (|l|), the difference between the positions (in gate voltage) of two consecutive 
peaks is 

6Vg{n) = {er])-'[U + - (2) 

Subtracting the constant charging term, the single-particle contribution to the energy spac- 
ing can be written as 

6Eg{n) = er] 6Vg{n) - U 

= En{Xn^i) — + 5En, (3) 

with SEn = En{Xn) — -E'„(X„„i). It is clear from Eq. (H) that the fluctuations in 5Eg{n) 
should not obey the Wigner-Dyson statistics when 5En also fluctuates on a scale larger or 
of the order of the mean single-particle level spacing A = (i?„(X„_i) — _E„_i(X„_i)). To 
understand this point, let us consider that a gate voltage variation 5Vg, sufficient to add 
one electron to the dot, also deforms the confining potential by 5X. Three illustrative 
situations are depicted in the single-particle spectrum shown in Fig. 1. The states indicated 
by the black dots fulfill Eq. (|ip for different values of n. The sequence (a) corresponds to 
the filling of single-particle states when 5X = 0, in which case the level spacings should 
indeed follow the Wigner-Dyson statistics. When 5X ^ and the deformation is moderate, 
we may have instead a sequence like (b). Now each filled state is slightly shifted with 
respect to its predecessor, level repulsion is weakened ,and deviations from the Wigner- 
Dyson statistics begin to appear. Finally, when deformations are strong, a sequence like (c) 
may occur: Consecutive filled states fluctuate independently and there is hardly any sign of 
level repulsion. In fact, in the presence of an overall downward drift, level spacings can be 
"negative" . 

The quantity which determines the form of the distribution of 5Eg, P{6Eg), is the ratio 
between 6X and the typical distance between anti-crossings in parameter space, X^. In what 
follows we study P{6Eg) function of SX/X^. 

Peak-to-peak fluctuations in 6X are of the order A/U and can be neglected when U ^ A, 
which is the regime relevant to the experiments. Thus, let us assume that 6X independs 
of n. This allows us to calculate P{5Eg) analytically in the perturbative regime. For this 
purpose, we model the parametric dependence of the electronic Hamiltonian as 

H,,,{X + 6X)=H + 6X K, (4) 

where H and K are large, independent N x N matrices belonging to the proper Gaussian 
ensemble. The orthogonal ensemble (GOE, /3 = 1) is used when no magnetic field B is 
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present, whereas for i? 7^ we use the unitary ensemble (GUE, j3 = 2). We use first 
order perturbation theory to describe the regime of small 5X/Xf. and write 
Therefore, 



5XK 



H,K' 



(5) 



where (■ ■ ■)h,k represents the ensemble averages over H and K. The average over K yields 
a Gaussian factor that depends on the first two moments of Knn- We choose {Knn) = and 
{{KnnY) = ^1^1- The former prevents level drift, while the latter corresponds to identify- 
ing the typical distance between anti-crossings with the inverse root mean square derivative 



of the energy levels with respect to deformations |T^, namely, X^. = A/ y {{dEn/dXy) . In 
these terms, 



P{6e, 



1 



27CX 



exp 



2^2 



(6) 



H 



where x = 6X/Xc. For the sake of simplicity, we have rescaled all energies in units of A, 
introducing 6eg = 6Eg/A and = En/ A for all n. Since Eq. (H) depends only on the 
spacing between neighboring levels, the average over H can be replaced by the convolution 



P{Seg 



dsP, 



WD I 



2tix 



exp 



2x2 



(7) 



P^^{s) denotes the Wigner-Dyson nearest-neighbor spacing distribution (WD), correspond- 
ing to the GOE (GUE) for /5 = 1 (/3 = 2) 0]. Notice that in the limit of a; ^ the Gaussian 
factor becomes a (5-function and P{Seg) equals the WD, as expected. 

The analytical evaluation of Eq. is possible when we approximate the WD by the 



Wigner surmise P, 



WD I 

/3 



the GOE, introducing = ^1/4 _)_ 1/(2x2), the result is 



with ai(2) = 7r/2(32/7r2) and 61(2) = 7r/4(4/7r). For 



P 



GOE 




exp 




whereas for the GUE, with ^2 



4x4^? 

= 4/7r + l/(2x^), we have 
2^/2fe, 



(8) 




7r5/2x3e| 

X exp 



(9) 



Notice that the Gaussian factor in Eq. (|^ has the unsatisfactory feature that it becomes 
indefinitely broad for increasing values of x. This is so because the perturbation theory 
cannot describe consecutive anti-crossings that stabilize the level position and therefore is 
limited to x <C 1. 



4 



We have also performed extensive numerical simulations of parametric realizations of 
Gaussian ensembles of = 500 random matrices to obtain P{Seg). For x < 0.25 the 
agreement with the analytical results is very good, with deviations becoming larger with 
increasing x, as expected. Moreover, only for x < 0.25 we have found that variations of x 
change significantly the shape of P{6eg). The numerical results indicate that, for x > 1, 
P{Seg) approaches a Gaussian distribution whose variance saturates slowly as x increases, 
as can be observed in Fig. 2. Thus, the utility of the analytical results is that they capture 
the essence of the change in shape of P{6eg), ceasing to be quantitative for large values of 

X. 

The reason why P{6eg) tends to a Gaussian may be explained in the following way. 
First, for a sufficiently strong deformation, we expect -E'„(X„) and £'„_i(X„_i) to become 
independent. Second, we have empirically found from our simulations that the position 
of each level taken alone fluctuates according to a Gaussian distribution, with a certain 
ensemble- dependent standard deviation (T/3 ~ 0(1). In view of Eq. (P|), we conclude that 



P{6eg 



1 



exp 



4a| 



(10) 



^47rcr^ 

(The variance of the unitary case is always smaller than the orthogonal one because the 



level repulsion grows with f3.) We remark that some analytical work [^, as well as the 
Dyson Brownian motion model hint the idea that the position of each level taken alone 
is Gaussian distributed, although we found no rigorous proof in the literature. 

So far we have not considered the possibility of having spin degenerate eigenstates. In 
other words, the discussion focused on the case that, due to some strong exchange contri- 
bution, the mean-field solution of the many-body problem led to states with broken spin 
degeneracy. Let us discuss now the opposite case of spin degenerate levels, when the nearest- 
neighbor level spacing distribution is 



pri^) = ls{s)+ipr{s). (11) 



For SX = this leads to a bimodal sequence of peaks spacings (and also of peak heights). 
However, in analogy to the previous discussion, as SX increases, even spin degenerate peaks 
will correspond to different dot shapes. Consequently, the bimodal structure will be rapidly 
destroyed. As 6X approaches the scale of one avoided crossing, P{6eg) will tend to a 
Gaussian. For x > 1 the only significant difference between the case of broken and preserved 
spin degeneracy is the variance of the distributions with respect to A, as shown in Fig. 2. 

We test our mechanism by modeling the dynamics of an electron in a quantum dot 
using the conformal billiard [^. The shape of the billiard is determined by the image 
of the circle of unit radius in the complex plane z under the conformal mapping w{z) = 
{z + bz'^ + ce'^z^)/{l + 26^ + Sc^)^/^. We studied the case b = 0.2 and c = 0.2, sweeping 
the parameter 6 in [7r/2,7r]. For this parameter range the conformal billiard is known to 



display chaotic motion in its classical limit ||21[. Another convenient feature of this model is 



the readiness to compute eigenvalues and eigenfunctions when the billiard is threaded by an 
Aharonov-Bohm flux of arbitrary strength. This allows us to study the cases of preserved and 
broken time- reversal symmetry with equal numerical effort. We have considered eigenvalues 
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ordered in ascending energy ranging from the 200**^ up to the 250**^. In this energy window, 
we followed the levels in parameter space to construct our histograms, observing that A and 
Xc stay almost constant over the whole energy-parameter window [^. The scaled results 
are in good agreement with our theoretical predictions based on RMT. In Fig. 3 we illustrate 
this statement showing the case of spin resolved degeneracy. 

Three experiments UToHr^ have obtained peak spacing distributions that fit Gaussian 
curves rather well. Our results are consistent with these findings if we take a; > 1, in- 
dicating that in the experimental setups there is at least one single-particle anti-crossing 
when the gate voltage is swept between peaks. Some difficulties arise when we try a more 
precise quantitative comparison. In particular, in Ref. [|1^] the authors found weak upward 
deviations from the Gaussian behavior at both tails of the distribution which cannot be 
reproduced within RMT and are not seen in our simulations. While the spectral statistics 
of the conformal billiard is rather universal, it is customary for chaotic systems to exhibit 
strong scarring and other non-universal features, resulting in levels that fluctuate weakly 
with shape deformations and do not follow RMT statistics. Although such levels are few, 
they may lead to deviations in the tails of P{5eg). This subject is currently under investi- 
gation. Thermal smearing was also not taken into account in our analysis and may lead to 
a small distortion of the peak spacing distribution. 

In conclusion, we have shown that the statistics of Coulomb blockade peak spacing 
distribution in chaotic quantum dots depends on shape deformations induced by variations 
in the gate voltage. Even when no fluctuations in the charging energy are present, provided 
that deformations are strong, the distribution of peak spacings is nearly Gaussian and quite 
distinct from the Wigner surmise, in qualitative accordance with experiments. We propose 
that the absence of spin- degeneracy bimodal structures seen in the experiments may also 
be caused by deformations. Although we have only considered the situation of constant 
charging energy in this letter, we believe that deformations should also be incorporated 
when wave function scarring and many-body effects beyond mean- field are present. 

We thank H. Bruus for providing us with the code for diagonalizing the conformal billiard. 
We gratefully acknowledge the financial support of FAPERJ, CNPq, and PRONEX (Brazil). 
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FIGURES 




gate voltage 
figure 1 

FIG. 1. Single-particle spectrum as a function of the gate voltage (or dot shape deformation). 
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figure 2 



FIG. 2. Standard deviation of P{S£g) as a function of x. Squares correspond to GOE and stars 
to GUE for the spin degenerate case. Triangles correspond to GOE and circles to GUE for the 
spin resolved case. 
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figure 3 



FIG. 3. Comparison between P{6eg) obtained from the conformal billiard (histograms) and our 
numerical simulations (solid lines) for the spin resolved case: (a) GOE for x = 0.25 (in gray) and 
x = 1 (in black); (b) the same for GUE. 
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